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CRITICAL MACH NUMBERS FOR 
SWEPT-BACK WINGS* 


by 


S. NEUMARK, Tech.Sc.D., A.F.R.Ae.S., A.F.I.Ae.S. 
(Aerodynamics Departmem, Royal Aircraft Establishment) 


INTRODUCTION AND SUMMARY 


The effect of yawing a wing in high-speed flow is to delay the onset of shock 
waves and to increase the critical Mach number. This is because shock waves can 
only develop along the isobars (running parallel to the edges of the wing, which is 
assumed untapered and infinitely long), and therefore only the velocity component 
normal to the edges is significant. To some extent the same is true for a finite sheared 
or swept-back wing, though the problem is made more complicated by the various 
additional factors, such as finite aspect ratio (tip effect), plan form (e.g. taper effect) 
and, what is perhaps the most important, the central kink effect. The experimental 
technique is extremely cumbersome because of the many shape parameters involved 
(see Figs. 1, 2, 3), and since analytical solutions present fundamental difficulties 
and pitfalls, designers tend to favour rule of thumb methods. The present paper 
contains a general review of the theoretical work done at The Royal Aircraft 
Establishment, which was limited to the case of zero incidence. This limitation is 
not a severe one, as flight at high speeds often implies low incidence. Incidences 
which are not negligible involve the solution of another fundamental problem (that 
of the lift distribution). 


The results presented here do not depend on lift theories, and belong entirely 
to the field of classical hydrodynamics. The method used is that of sources and 
sinks and is really a linear perturbation theory for thin aerofoils. The early results 
obtained were confined to the kink effect and referred primarily to infinite aspect 
ratio. It is felt now that the fundamental mathematical difficulties connected with 
the kink, tips and taper have been satisfactorily overcome, and the results are shown 
in the series of curves presented. These give velocity (and hence pressure) 
distributions, the course of the isobars, and critical Mach numbers as functions of 
geometric parameters. For a finite swept-back wing the central kink section is 
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Fig. 3. 
Tapered swept-back wing. 
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CRITICAL MACH NUMBERS FOR SWEPT-BACK WINGS 


normally the first danger zone, and the lower critical Mach number is defined as 
that corresponding to the flow reaching sonic characteristics locally in this section. 
With increasing Mach number the sonic region gradually expands spanwise, and 
the upper critical Mach number is the one above which shock waves may cover 
almost the entire wing. The tips are often immune up to even higher speeds. 


The effect of sweepback is of great practical interest since it seems to make 
possible flight at very high subsonic or even low supersonic speeds, while retaining 
essentially subsonic aerodynamic characteristics. 


Notation 


F (x) 

F’ (x) 
f(x, y) 
f’.(x, y) 
G 


aspect ratio of complete swept-back wing (or double the aspect 
ratio of a finite sheared wing) 


speed of sound 

speed of sound in undisturbed flow 
half-chord of wing section (Fig. 1) 
local drag coefficient in potential flow 
chord of wing section 

root chord (Figs. 2 and 3) 

tip chord (Figs. 2 and 3) 

function determining wing section 
first derivative of F (x) 

function determining wing surface 
first partial derivative of f(x, y) with respect to x 


ratio of the true local supervelocity at any point of the surface 
of an untapered swept-back wing to the maximum supervelocity 
on an infinite sheared wing of the same section and angle of 
sweepback 


value of G for the “ Géthert wing” 

maximum value of G 

Mach number of undisturbed flow (or flight Mach number) 
critical Mach number 

local source strength 

auxiliary symbols, defined in (17) 


semi-span 
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thickness of wing section 
value of t for the “ Géthert wing ” (Fig. 23) 
undisturbed velocity of air flow (velocity of flight) 
resultant velocity on wing surface 
velocity component normal to the leading and trailing edges of 
an untapered sheared or swept-back wing 
velocity components induced by a system of sources and sinks 
representing a wing 
rectangular co-ordinates (x chordwise, positive forwards; y span. 
wise, positive to starboard; z vertical, positive upwards) 
chordwise co-ordinate of source filament 
x*, y* co-ordinates of a source element in the x, y-plane 
Y=y/s auxiliary variable (see (24)) 
y adiabatic constant (~ 1.4) 
8=|v,/U |max | maximum supervelocity ratio in incompressible flow 
value of 6 for the “ Géthert wing” (Fig. 23) 
coefficient of convergence of a tapered wing (Figs. 2 and 3) 
ratio of the section area to that of a circumscribed rectangle 
non-dimensional co-ordinates in the x, y-plane (see (17)) 
value of € in the central kink section at the point where the 
supervelocity attains its maximum 
angle of sweepback 
angle of sweepback of mid-chord line (Fig. 3) 
angle of sweepback of leading edge (Fig. 3) 
angle of sweepback of trailing edge (Fig. 3) 
angle of sweepback of the “Géthert wing” (Fig. 23) 
taper ratio 
thickness ratio of wing section 


PART I 


VELOCITY DISTRIBUTIONS AT ZERO INCIDENCE IN 
INCOMPRESSIBLE FLOW 


The problem of finding the velocity distribution on a swept-back wing being | 
essentially three-dimensional, there is little hope of solving it exactly by potential 
flow theory. The only suitable approach seems to be the approximate first order 
method (usually called “linear perturbation theory”) in which the wing is replaced 
by an approximate continuous system of sources and sinks, and the velocities [ 
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CRITICAL MACH NUMBERS FOR SWEPT-BACK WINGS 


induced by such a system are calculated. The field of these induced velocities must 
then be superimposed on that of the undisturbed uniform velocity U. If the induced 
velocity components are denoted by v,, v,, v., the resultant velocity V is given by 


(1) 


+2,? +2,’. 


All three induced velocity components are small of the first order, and the resultant 
velocity, to the same order of accuracy, is 


V ~-U+2,. (2) 


This is the value to use when finding either local pressure or the shape of the 
isobars, but for determining critical Mach number the component v, may be as 
important as vz. 


The linear perturbation theory has been applied successfully by many authors 
to the two-dimensional case, i.e. to a straight wing of infinite aspect ratio. In this 
case the local source strength g (per unit area) at any point of the chord is 
proportion to the local slope of the wing section, that is 


q= —2UF’ (x). (3) 


where z=F (x) is the equation of the upper arc of the wing section. 


The infinitesimal induced velocity components dv, and dv, contributed by one 
single source filament are shown in Fig. 4. The total induced velocity components 
are obtained by integration, and it can be proved that for a thin wing, the total v,, 
whether calculated at a point P on the surface or at its projection Q on the chord, 
is, to the required order of accuracy, 


lim v,= —UF’(x). (4) 
230 


Hence the boundary conditions are satisfied, to the same order or accuracy. The 
supervelocity v, is then determined by integration, and it is again sufficient to 
calculate it at the point Q using the formula 


b 
F’ 
x— 

b 


IS 


The integral involved is improper, and its principal value should be taken. Thus the 
supervelocity distribution is determined, in the most general case, directly from the 
equation of the section. 
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The example in Fig. 4 refers to the simplest “bi-convex parabolic section,” 
given by 
2 
z=F 


The source and sink distribution is linear in this case (Fig. 4), and the supervelocity 
becomes 


gt 


where #=t/c is the thickness ratio, and £=x/b is the non-dimensional chordwise 
co-ordinate, which is equal to —1 at the trailing edge and +1 at the leading edge, 
The chordwise supervelocity distribution is shown in the lower diagram of Fig. 4. 


Similar solutions can be obtained for an arbitrary section shape. Their accuracy 
is quite good for thin sections, though singularities usually occur at both the leading 
and trailing edges. The important maximum values of v, obtained by this method Ff 
are sufficiently accurate for all practical purposes. In the given case | v.|max Occurs 
exactly at the mid-point of the chord, and the maximum supervelocity ratio is 


= 4 1.278. (8) 


max 


For an elliptic section it may be noticed that 5=#, this being almost the minimum 
supervelocity ratio obtainable, and, from the point of view of critical Mach number, 
it is important to keep the ratio 5/? as small as possible. It will always exceed 
unity and will usually be greater than 1.2. 


If 5 is determined for any section, the maximum total velocity is 
V max = U (1 +8), (9) 


and the condition for critical Mach number is obtained by equating this expression 
(transformed to allow for compressibility) to the local speed of sound. (For the 
results see Part II). 


Since the method of sources and sinks is so effective in the two-dimensional 
problem, it is logical to try to apply it to three-dimensional cases, especially for 
swept-back wings. In the general case of a thin wing with arbitrary plan form and 
arbitrary distribution of (symmetrical) sections, the plan form area S will be defined 
in terms of x and y, while the upper surface of the wing is represented by the equation 


z=f (x, y), . (10) 
(x, y within S). 
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FILAMENT 
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DISTRIBUTION 
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4-0-5 


Fig. 4. 
Distribution of sources and induced velocity for a straight wing with 
bi-convex parabolic section. 
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It may then be assumed that the area S is covered by a continuous system of source; 
and sinks, their local intensity per unit area being 


(x*, y*, are the co-ordinates of an arbitrary source element). The v, component 
of the induced velocity at any point (x, y) is 


(x -— x°). f° (x*, y*) 
2 (x—x* +(y—y*)? 


dx* dy*, . (12) 


and similarly for vy. It can be proved that this generalised method still gives 
solutions accurate to the first order of thickness ratio, although the errors may be 
somewhat larger numerically than in simple two-dimensional problems. The method 
fails only at the contour of the plan form, which is a locus of singularities on which 
either v, or v, or both become logarithmically infinite. This is not serious, as it 


affects only very narrow marginal strips, and the velocity distribution over almost | 


the entire wing surface is represented with sufficient accuracy for practical purposes. 


For simplicity the bi-convex parabolic section shown in Fig. 4 has been used 
in the first instance, because it gives a linear chordwise source distribution and is 
therefore a perfect “guinea pig” section. The calculations can be made for any 
arbitrary section. 


The simplest case to which the generalised method has been applied is that of 
an infinite sheared wing (Fig. 21). This is an interesting test case because the 
solution may be obtained either by a direct application of the general method of 
sources or by resolving the undisturbed uniform flow into components along and 
normal to the wing edges, thus reducing the problem to a two-dimensional one. 
Both methods give the same results, i.e. 


4 l 
v.=— 108.5) cos ¢ 


€ is still the convenient non-dimensional chordwise co-ordinate, becoming —1 at 
the trailing edge and +1 at the leading edge. In this case the isobars are all straight 
lines parallel to the edges, and the maximum relevant velocity, needed for deter- 
mining the critical Mach number, is the maximum normal component (at mid-chord 
point in the present example) 


(Vi)max=U (cos¢+8). (14 
92 
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Comparison of this relation with equation (9) shows that, because of the presence of 
the term cos 9, (Vn)max can be reduced by increasing the sweepback. 


The next step in the theory is to apply the general method to the case of a 
swept-back wing of infinite aspect ratio (see Fig. 1 with s=00), as first suggested by 
H. Ludwieg'”) of Géttingen. His method was essentially the same as that of the 
present paper, but unfortunately it led him to a false result. He assumed rightly that 
for such a wing the danger zone was at the kink, and he confined his calculation to 
this section. A simple integration brought him back to formula (13), exactly as for 
an infinite sheared wing, but he failed to notice the full paradoxical meaning of 
his result, and used it to calculate the critical Mach number in the danger section. 
According to Ludwieg, the maximum total velocity in the kink section would be 


V | max = U (1+8 cos ¢), ° (15) 


8 being the maximum supervelocity ratio in the two-dimensional case. 


Ludwieg’s error remained unnoticed until calculations of the velocity 
distribution over the whole wing area were made at the Royal Aircraft Establish- 
ment. The result was the following, regrettably complicated, formula 


1+£&+(r,—n)tane 
1—£—(r,—n)tan¢ 
r,+(1+€&)sin ¢ cos ¢ +n cos 
r,—(1—)sin¢ cos ¢ cos 2) 


4 U9 cos ¢ {1 log. 


+4n tan ¢. log. 


(17) 


The corresponding formula for v, was still more complicated. 


Assuming »—>0 we obtain the true formula for the velocity distribution at 
the kink section 


1—€ 


This differs from (13) by the factor (1 +sin ¢)/(1—sin¢) inside the bracket. 


The difference is due to Ludwieg’s assumption that »=0 before the final 
integration, while the reverse procedure is adopted here. A mathematical singularity, 
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S. NEUMARK 


connected with the geometrical discontinuity in plan form, implies that “ the limit 
of an integral” is not equal to “the integral of a limit.” Because of this Ludwieg 
missed what we now call “kink effect.” 


The final v, distribution over the whole wing is shown in Fig. 5. It may be 
noticed how this varies continuously spanwise from the symmetrical picture at 
n=00 to a completely asymmetrical one at »=0 (kink section). The maximum 
supervelocity is located farther back, and its numerical value is much higher than 
that obtained by Ludwieg. By introducing the function 


G()= 1-4 log. 


1—€ 1-—sing (19) 


the total velocity at the kink section may be represented by 
| V |=U (1+8G cos ¢), 
and its maximum by 
| V | max =U (1 + 8Gimax COS 9). 


The factor Guax, which depends on ¢, shows the difference between Ludwieg’s 
result and the present one. 


The isobars on the swept-back wing of infinite aspect ratio are shown in Fig. 6. 
The kink effect can be seen clearly, and at some distance from the kink the isobars 
become almost parallel to the edges, the conditions gradually becoming identical 
with those on an infinite sheared wing. 


Figure 7 shows the function G for different angles of sweepback, and Fig. 8 the 
product Gcosy. The kink effect increases rapidly with ¢, but the final effect of 
sweepback is favourable, though less than that estimated by Ludwieg. Fig. 9 shows 
this more clearly, the curve cos ¢ corresponding to Ludwieg’s estimate. G,,,, shows 
the error, and Gx cos ¢ the true effect of sweepback in the danger section. 


The results above are limited in their application as they are confined to the 
bi-convex parabolic profile, which is never used in practice, but they can quite 
easily be generalised for arbitrary profiles, particularly as regards the kink section. 
The general formula is as follows 


= { (x) log. COs ¢, (22) 


where F”’ (x) denotes the first derivative of the equation of the section. The velocity 
distribution at the kink may therefore be determined immediately from that of the 
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given section in two-dimensional flow. An example is given in Fig. 10 for a 
round-nosed profile. 


The series of figures which follow show the results of similar calculations for 
various plan forms. Figs. 11 and 12 show the velocity distribution and the isobars 
on an untapered swept-back wing of finite aspect ratio. The kink effect is seen to 
be almost exactly the same as in the previous case, but in this case a new 
phenomenon, the tip effect, makes its appearance. There is a strong tendency for 
the supervelocities to decrease towards the tips, while the position of the local 
maximum moves forward and the isobars bend near the tips to become almost 
parallel to the flow. This seems to indicate that at zero incidence the tips are the 
least prone to shock-stalling. A warning, however, is now necessary. As mentioned 
previously, the accuracy of the method is not adequate within a narrow marginal 
strip, and therefore complete agreement between theory and experiment is hardly 
to be expected very near to the tips. This tip effect may well prove troublesome in 
the future, especially if the tips are square-cut or provided with end-plates. With 
large end-plates, the tip velocity distribution would be quite similar to that at the 
kink, but reversed fore and aft. For well-rounded tips the results presented here 
should be a good approximation, although really accurate results can only be 
obtained by more refined methods. 


The interesting point is that, even with the small aspect ratio assumed (A =2), 


the isobars run almost parallel to the edges over a large part of the wing and the 
velocity field is almost the same as on an infinite sheared wing. For determining 
critical Mach numbers for wings of fairly large aspect ratio, only two types of 
section need be considered : 


(i) The kink section, which is the first danger zone. 


(ii) The “quasi-two-dimensional” section, i.e. a section representative of the 
large area where conditions are similar to those on an infinite sheared 
wing, and may therefore be examined on the basis of two-dimensional 
theory. When critical conditions appear here, the entire wing may be 
considered as shock-stalled. This point is considered in the second part 
of this paper. 


Figures 13 and 14 show the velocity distribution and the isobars on a 
semi-infinite sheared wing, which is a particularly simple case for studying the 
isolated tip effect. Figs. 15 and 16 show similar pictures for a finite sheared wing. 
One point in connection with these last figures may be noted. The maximum 
supervelocity in the centre of a finite sheared wing might be expected to be slightly 
smaller than on an infinite sheared wing, due to the cutting off of the infinite “ tails,” 
carrying their systems of sources and sinks. However, it can be shown that, if 
sing >1//3, or ¢>35° 16’, the effect is opposite because of the peculiar 
disposition of sources and sinks. This is shown by the number 1.047 in Fig. 16, 
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UNSWEPT WING 


CENTRAL SECTION, 
SWEPT-BACK WING 


Round-nosed profile, sources and supervelocity distribution. 


which means that the maximum supervelocity exceeds by 4.7 per cent. that obtained 
for an infinite sheared wing. This surprising effect is shown in more detail in Fig. 17. 


There is a very similar effect of the finite aspect ratio of swept-back wings on 
the maximum supervelocity in the kink section. An approximate formula, which is 
sufficiently accurate if the aspect ratio A is not too small, is 


(3 sin? ¢ — 1) cos* ¢, ‘ - (23) 


where v denotes the ratio of the profile area to that of a circumscribed rectangle. 
The effect may usually be neglected. 


In the case of a tapered wing, there are further serious difficulties. First, the 
integrals involved become increasingly complicated and can no longer be reduced 
to elementary or tabulated types. Second, the number of parameters involved 
increases alarmingly (angle of sweepback ¢,, coefficient of convergence «, and taper 
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Fig. 17. 
Maximum supervelocity in the central section of a 
finite sheared wing. 


ratio yf). In spite of these difficulties the method can still be used, but the 
computation involved is considerable. Only the very simplest case of an unswept 
tapered wing, especially that of zero taper ratio, ie. a rhombus wing (Fig. 2) is 
considered here. The maximum supervelocity at the centre can be represented by 
the formula 


log. (1 - dY (24) 


1 
4 


(Y=y/s). 


More complicated formulae have been obtained for the more general cases. The 
results are illustrated in Fig. 18, which shows the effects of the parameters « and Y. 


Or, alternatively, ~,, ¥ and A (aspect ratio); the parameters used in the text are more convenient. 
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i 
Fig. 18. 
Variation of maximum supervelocity of straight tapered 
wing (at mid-point of central section) with taper ratio y 
and coefficient of convergence ¢; curves of constant aspect 
ratio A shown as broken lines. 


The effect of « is surprisingly large, but the effect of ¥ is comparatively insignificant 
until Y becomes quite large. However, large values of ¥ correspond to absurdly 
small aspect ratios, as shown in the diagram. 


Similar results for the most general plan form of Fig. 3 are not yet complete, 
and to represent them adequately a whole series of diagrams would be needed. 
However, the results presented give a fairly good idea of the significance of the 
various shape parameters. 
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One further point is of interest. The drag of any finite body in potential flow 
is always zero, but for an infinite swept-back wing or a semi-infinite sheared wing, 
the total drag has finite values, which increase with ¢ up to a certain maximum, 
and then fall again. This is due to the kink effect or to the effect of an oblique 
tip, both of which destroy the fore-and-aft symmetry of the velocity distribution. 
Fig. 19 shows the drag distribution for both types of wing, and it is remarkable 
that the total drag (area under the curve) is the same in each case. In this connection 
it may be remarked that the total drag of a finite swept-back wing is zero, although 
local drag exists in the kink region and local thrust in the tip regions. 


PART II 
DETERMINATION OF CRITICAL MACH NUMBERS 


Part I of this paper dealt with the determination of maximum velocities, either 
for the entire wing or for particular sections, in potential incompressible flow. 
The purpose of Part II is the evaluation of critical Mach numbers, and the results 
of Part I are used to find the corresponding maximum velocities in compressible 
flow. There are several ways of doing this, but only the simplest Glauert-Prandtl 
method, with Géthert’s generalisation for three-dimensional problems, is used here. 
More refined and more recent methods can be applied (e.g. the Karman-Tsien 
method) but they require much more computational work. 


The first step is to determine the critical Mach numbers for any given wing 
profile in simple two-dimensional flow. The fundamental condition is that the 
maximum velocity just equals the local speed of sound. V,,.; in incompressible 
flow is given by (9); in compressible flow, it becomes 


=U (14+ ‘ ‘ & 


where M=U/a, is the Mach number and a, is the speed of sound in undisturbed 
flow. The local speed of sound at the critical point can be determined from the 
generalised Bernouilli equation as 


2 1-—M? 


Equating (25) and (26), gives the fundamental relationship 


between maximum supervelocity ratio and critical Mach number in two-dimensional 
flow. This relationship is represented in Fig. 20 over a large range of 5. The range 
of practical importance is between 5=0.05 and 0.2, and in this interval the derivative 
dM / dé differs little from unity. A rough mnemotechnic rule follows that a reduction 
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of 0.01U in maximum supervelocity gives an increment of about 0.01 in critical 
Mach number. 


The next step is to determine M., for an infinite sheared wing. It is important 
to differentiate between a sheared and a yawed wing, and the differences are shown 
in Fig. 21. Most experimental results refer to yawed wings because of the obvious 
ease of rotating a model in a wind tunnel through a suitable angle, instead of making 
separate models for different angles of sweepback. This has often led to confusion 
when comparing theoretical and experimental results. 


The methods of treating each case are very similar and consist of considering 
only the component flow normal to the wing leading and trailing edges, while the 
Mach number refers to the full undisturbed velocity U. The local speed of sound 
at the critical point is now 


? . . 2 


and this must be equated to the maximum normal velocity component on the wing 
(cf. the formulae in Fig. 21). The relationships between 4 and critical Mach 
number follow. 

1 — M? cos? ¢ | 


For yawed wings: 5=4/(1—M?* cos? ¢) “Mad; 


(29) 


M’ cos? ¢ 


for sheared wings: 8= /(1—M?*cos?¢) 1} cos 


(30) 


The two formulae are illustrated in Fig. 22 which shows that the differences 
between the two cases are considerable and increase with increasing 5. The most 
important point to note is that with large angles of sweepback, critical Mach 
numbers can be obtained which are well above unity. If it were possible to design 
wings possessing such properties over their entire surface, aircraft would be able to 
fly at supersonic speeds while preserving subsonic aerodynamic characteristics. This 
was suggested by Busemann at the 1935 Volta Congress. 


The weak spot of a true swept-back wing (Figs. 1 and 3) is the kink section and 
the surrounding region. Here the isobars must always be at right angles to the 
direction of flight, the relevant maximum velocity (Fig. 23) is always greater than U, 
and therefore the critical Mach number at the kink section must always be less than 
unity. This may be called “lower critical,” while that corresponding to the main 
“ quasi-two-dimensional” region will be the “upper critical” Mach number. 


The velocity field in compressible flow for a complete swept-back wing can be 
determined by considering a hypothetical wing (“ Géthert’s wing”) in incompressible 
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SHEARED WING 


YAWED WING / 


Fig. 21. 


Velocity components on infinite sheared or yawed wing, relevant for determining 
critical Mach number. 
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SHEARED WING 
YAWED WING 


UPPER CRITICAL 
———~ LOWER CRITICAL 


Critical Mach numbers for infinite sheared Lower and upper critical Mach numbers 

and yawed wings of varying angle » against of swept-back wings with bi-convex 

maximum supervelocity ratio in two-dimen- parabolic section for varying sweepback 
sional flow 38. angle ¢ and thickness ratio #. 


flow, as shown in Fig. 23. This method is especially useful for the kink section, 
while the conditions in the “quasi-two-dimensional” region may be examined 
either by Gothert’s method or simply by resolving the flow as for an infinite sheared 
wing. Both methods have been used and lead to identical results. 


Assuming that the aspect ratio is not small, the maximum velocity at the 
kink section is 


8G’ max COS 
(1 — M? cos? ¢) 


| (3) 


where G’max (the kink effect factor, see Fig. 9) corresponds expressly to the angle 
of Gothert’s wing. 


Equating (31) to the relevant local velocity of sound, still given by (26), we 
obtain for the lower critical Mach number 
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This formula can be used only with reference to the particular wing section, because 
the factor G’n.x depends on the section shape as well as on g and M. Fig. 24 shows 
the results for wings with bi-convex parabolic section, and both the upper and 
lower critical Mach numbers are shown as functions of sweepback angle ¢ and 
thickness ratio #. For large ¢, there is a wide gap between the two critical values. 
When the true Mach number increases from the lower to the upper critical, the 
supersonic region gradually spreads out from the centre, eventually covering the 
entire wing. The calculation of the effect of this process on the wave drag is beyond 
the scope of this paper. 


The lower critical Mach number for tapered swept-back wings (particularly 
delta wings) can be calculated in a similar way, but this requires a great deal of 
computational work. There is difficulty in defining the upper critical Mach number, 
since the simple method of resolving the flow into two components cannot be 
applied directly. The question of what should be considered as the mean or 
representative angle of sweepback is still an open one. 


There are many more problems to be solved in this field, the most important 
being to extend the theory to cover wings at variable incidence. 
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GENERAL PERFORMANCE REDUCTION 
EQUATIONS FOR RECIPROCATING-ENGINED 
AIRCRAFT WITH CONSTANT-SPEED PROPELLERS 
by 


W. J. D. ANNAND, B.Sc., A.F.R.Ae.S. 
(Rolls-Royce Limited, Hucknall) 


SUMMARY 

The method of analysis adopted by Cameron"? has been extended to develop 
more general equations for the correction of level speed and rate of climb for 
variation of air temperature and aircraft weight, applicable to aircraft with 
reciprocating engines of any type driving constant-speed propellers, operating at 
speeds sufficiently high for compressibility effects on drag and propeller efficiency 
to be appreciable. Allowance is made in the equations for ram effect on 
engine power. 


The practical application of the equations is discussed, and illustrated by a 


worked example. The influence of errors in the estimation of the parameters 
expressing the variation of engine power and propeller efficiency on the corrections, 
in this particular case, is discussed. 


1. INTRODUCTION 

In order that measurements of the performance of an aircraft made at different 
times may be correlated and compared, it is necessary to have a method of 
determining or eliminating the effects of variation of those parameters which cannot 
be held constant. This process is called “performance reduction.” In what are 
known as “analytical” methods of performance reduction, the variation of 
performance with the uncontrolled parameters is determined explicitly, and the 
observed performance is reduced to standard values of the parameters by the 
numerical correction of individual observations. Inj practice, the uncontrolled 
parameters are most commonly the air temperature and the weight of the aircraft. 


Analytical methods of performance reduction for piston-engined aircraft with 
constant-speed propellers, operating at Mach numbers below the critical value for 
drag rise, were developed in 1942 by Cameron"), The working methods set out in 
that reference are in general use at the present time. Their application, however, 
is limited by the purposes for which they were developed: they were intended to 
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yield approximately correct results for all types of aircraft falling within the above 
category, and for this purpose certain simplifications were adopted and average 
numerical values of certain parameters were employed. 


In the present paper, the author has followed the general line of approach 
adopted by Cameron, to develop more general equations for the correction of level 
speed and rate of climb for variations of air temperature and of aircraft weight, at 
constant air pressure and propeller r.p.m. The equations given include allowances 
for the effect of compressibility on aircraft drag and on propeller efficiency, and for 
the effect of ram on engine power. They are developed in a form such that they 
may be applied to aircraft powered by any type of reciprocating engine driving a 
constant-speed propeller. 


As neither Cameron’s paper, nor any other on analytical methods of 
performance reduction, has yet been published, to the author’s knowledge, the 
derivation of the equations is given in detail, in the belief that the indication of the 
methods of attack applicable to this and similar problems will be of value to students. 


Notation 

A aspect ratio of wing 

b span of wing 
engine supercharger constant 
total drag coefficient 
Cy — KC,?/("A) 
value of C,, at low Mach number 
lift coefficient 
propeller power coefficient 
engine cylinder compression ratio 
total drag of aircraft 
propeller diameter 
propeller parameter =(J/ 1) 0n/0J 
propeller parameter =(M,/1) 0n/0M, 
propeller parameter =(C,/ 1) 0n/0C, 
P,-P 


[m* / [2m*+1] 


propeller advance coefficient = V;/(nd/o) 


(xA) times the slope of the curve of C, against C,? at constant 
Mach number 
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Mach number 

propeller tip Mach number 

Vi/Vima 

propeller rotational speed 

power delivered to propeller shaft 

indicated power of engine (i.e. power developed in cylinder) 
ambient air pressure 

stagnation pressure of air in engine intake 
wing area 

stagnation temperature of air in engine intake 
ambient air temperature 


standard value of air temperature to which the performance is to be 
corrected 


equivalent air speed =(true air speed) x /o 


equivalent air speed for minimum drag at low altitude 
=( 4K W? 0.25 


true (tape-measure) rate of climb 


value of v. corrected to the I.C.A.N. standard 
temperature at the test height 


rate of climb indicated by an altimeter calibrated according to the 
L.C.A.N. convention 


weight of aircraft 
4p 
K/(4%p,b’) 


engine parameter, depending on the rotational speed, representing the 
charge temperature rise from intake to cylinder 


deviation of temperature from standard =t — f, 
deviation of weight from standard 


correction to be added to observed V; to obtain the value at the 
standard temperature or weight 


correction to be added to observed v, to obtain the value at the 
standard temperature or weight 


Coz /Cozo 
M/M., 


ratio of the specific heat of air at constant pressure to that at constant 
volume 
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propeller efficiency 

air intake efficiency 

(M/Cyz)dCy,/dM 

(t/P)oP/dt 

(V;/P)aP/aV; 

air density 

standard value of air density at sea level 
P| Po 

211, (1+ Ep) +E, 


ASSUMPTIONS MADE 
The following basic assumptions are made: 
(i) The drag coefficient Cy of the aircraft may be expressed as the sum of two 


parts, one of which, Cpz, is independent of C,, and the other of which is proportional 
to C,?. In the usual notation, 


K 
Cus + C,?. 


This assumption is that generally employed in performance calculations. It is 
known to be inexact at high C,, but is considered to be sufficiently accurate for 
the calculation of small performance changes, throughout the normal flight 
range of speeds. 


(ii) Only the part Cp; of the total drag varies with Mach number. Evidence 
as to the variation of K is scanty and, for the present, the additional complication 
introduced by allowing for such variation does not appear justified. It is convenient 
to write Cpz=2 Cpz, where Cpz, is the value of Cpz at low (nominally zero) Mach 
number. 

(iii) The power P delivered to the propeller shaft is a function only of the air 
pressure p, the air temperature ft, and the equivalent air speed V;. The inclusion 
of V; in the parameters allows for the effect on the engine power of the ram pressure 
and temperature rises in the air intake of the engine. 

(iv) The propeller efficiency 4 is a function only of the advance ratio J, the 
power coefficient C,, and the tip Mach number M,. 

(v) No thrust is obtained from the exhaust gases. This assumption is made to 
avoid excessive complication in the analysis. A method of estimating the effect of 


variation of exhaust thrust, and some indication of the magnitude of this effect, are 
given in Ref. 2. 
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3. BASIC EQUATIONS 


Using assumptions (i), (ii) and (v), the well-known performance equation is 
obtained 


To simplify manipulation this is written in the form 


. (1) 


where Y and Z depend only on the aircraft configuration. The total drag of the 


aircraft is given by 


whence is obtained the equivalent air speed at which the drag would be a minimum 
if no compressibility effects existed, 


vO 
al 


V ima = (ZW? 


( 4KW? 
tp 


At constant air pressure and engine r.p.m. the following relations hold: 
ot OJ/\J Ot n 0C,/\C, Ot OM or 


E, 1 oP 


Pot co Of 


= (Es + (I +II,)—(1— 8) Ey] 


Vv, [E, + Epil, + 


Md 


is 
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4. VARIATION OF PERFORMANCE WITH TEMPERATURE AND 
WEIGHT AT CONSTANT PRESSURE AND ENGINE R.P.M. 


4.1. Vanation of sheed in level fught with temperature 


Putting v.=0 in equation (1) and differentiating with respect to temperature 
at constant weight, 


24V, _ ZW* dV, 


ndt Jo dt YoV3+ZW?2/V, 
Substituting from the relations given in Section 3, 


ld 
[Es + + BEw] + + (1+ (i B)Ey) + 


1 dv, 


V4 1] dt 


am* +1 Vi dt’ 


which gives on re-arrangement, 


t dV, 4201, 
fit f (m, M)—(E, + BEy)— Uy (1 + Ep) 


4.2. Vanation of speed in level fight with weight 


Similarly, by differentiation at constant temperature, 


dz dV, ,dV, ZW? dV, 
YV; dV, +3YaV, dw aw + 


which gives, on substitution, 


(zu +32)-1 


a+] 


i 


| 
: 
V, dt 2t Bi 
OV, dW” P OV, dW” 
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and on re-arrangement 


V, dW  (am* + 1)[f (m, M)—(E, + BEu)— (1+ Ep)] 


4.3. Variation of rate of chmb with temperature at constant V, 
Differentiating equation (1) at constant Vi, 


Ot P Ot vo dt 


which gives on substitution 


% 
1 dt 2t 


Jo ZW YV; 


and therefore 


V; RC 


1 2 2 


x [2m (1 +Ep)+E,; +2Ep (1 1| 


4.4. Variation of rate of cimb with weight at constant V, 


In this case, differentiation with respect to weight gives 


dv./dW 


Therefore dv,|dW = -(1/W) [ 0) | 


5 WORKING EQUATIONS 

From the equations already developed the following working equations can be 
obtained for the corrections to be added to the observed values of V; and », for 
deviations At and AW from the standard values of temperature and weight. 


Level flight 


4oAt/t 
V; M)-(E, + BEy)— (1+ E,) 


AV, _ 2AW/W 
V, (am*+1)[f(m, M)-(E; + BEw)— Ty (1 + Ep)) 
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4K 
WU xp, 


where ¢ represents the sum 
211, (1+ Ep)+E,+2Ep—(1—8) Ey—1. 


It should be noted that ¢ will in general be negative so that AV, is positive for 
positive Ar=1t—1,. Note also that in the equations as presented here consistent units 
are employed. 


These equations reduce to those developed by Cameron”) for 2z=1.0, 
p=Ep=E, =II1,=0. 


The correction Av, given for a temperature deviation Aft is the correction to 
the true tape-measure rate of climb. Usually, the rate of climb is measured by the 
variation of the reading of an I.C.A.N.-calibrated altimeter. If a test is made ona 
day when the temperature is other than the I.C.A.N. standard, the altimeter does 
not indicate the true rate of climb but gives the variation v., of the pressure height 
with time. The relation between v,, and 2, is 


© 


where ¢ and ¢, are the test and standard temperatures at the pressure height in 
question. Thus 


and the correction to be added to v,, to give the true rate of climb v,, at the standard 
temperature f, can be written 


A V 0s 
Vea Veg [ (am? + 1/m*)] } 


neglecting terms of second order in At/t. 


Note that if it is desired to correct the observed rate of climb to conditions other 
than the LC.A.N. standard, this last equation cannot be used, and equation (4) must 
be employed, first correcting all observations to the true rates of climb by (6). 
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6. PRACTICAL APPLICATION 


The dependence of the “compressibility” parameters 2, « and Ey and of the 
ratio 8 on the design of the aircraft is such that it is impossible to make use of 
average values to construct simple general curves from which the corrections might 
be read off for any aircraft. However, for any given case it is a fairly simple task 
to construct sets of curves from which the corrections may be read off as functions 
of equivalent air speed and altitude, on the assumption that the substitution of 1, 
for t in the correction formulae will not seriously affect the accuracy, so that altitude 
may be taken to determine both the pressure (which with V; determines M) and the 
value to be used for ¢ in the formulae. This assumption is clearly valid for all cases 
in which the use of a “ differential” correction is justified. 


The construction of such curves will be illustrated by a numerical example 
(Section 7), but some general remarks on the derivation of the various parameters 
may be helpful. The various drag parameters Cp,., K, 2 and » may be obtained 
from wind tunnel measurements, or by analysis of flight performance measurements. 
a and » will be functions of M only. Experience suggests that in general the values 
of E, and Ep will vary little with flight condition within a given régime—i.e. climb 
or level flight— and that for a given régime they may be considered to vary linearly 
with V; and altitude; thus they need only be calculated for a few limiting conditions. 
The method of estimation of propeller efficiency recently published by the S.B.A.C.®) 
is very suitable for the estimation of all three propeller functions. The values of 
II, and II, will usually be fairly easy to obtain, either graphically from test-bed 
measurements of the response of the engine to variation of intake pressure and 
temperature, or analytically from formulae for the engine performance when these 
are available. If assumption (iii) holds, 11, and II, also are (with sufficient accuracy) 
functions of altitude and V; only. The derivation of formulae for Il, and II, for 
the conventional four-stroke mechanically supercharged engine is given in Appendix I 
as an illustration. 


7. ILLUSTRATIVE EXAMPLE 


As an example of the construction of working charts for a particular case, the 
corrections to level speed for weight and temperature variations from standard will 
now be worked out for a hypothetical aircraft, and the influence of certain of the 
parameters in this case will be discussed. The basic aircraft data are 


Wing area S=200 ft.? 
Wing span b=35 ft. 
Cp, =0.015 
K=1.20. 
The weight chosen as standard is 10,100 Ib., at which Ving is 200 m.p.h. 
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Fig. 1. 
Variation of a and » with M for example aircraft. 

In Fig. 1, 2 and m are plotted against Mach number. The Mach number at 
which the drag rise begins has been chosen somewhat lower than would be expected 
on a well-designed aircraft, in order to illustrate the effect of the drag rise on the 
corrections. 


The maximum Mach number for level flight at sea level has been taken as 
0°50. The power required for this performance is 3,500 horse power and the 
propeller diameter has been chosen to absorb this power according to the recom- 
mendations of Ref. 4. The propeller details are 

Diameter d= 16 ft. 
Number of blades=4 
Solidity at 0.7 radius=0.15 


Thickness /chord ratio=0.07 at 0.7 radius 
0.20 at root. 


The engine parameters are 
C=0.002 
z=150°. 
For simplicity, F/P has been taken to be constant at 0.10, and the intake efficiency 


y has been taken to be 0.90. The maximum engine speed is 1,194 r.p.m. at 
the propeller shaft. 


From these data, the values of the engine and propeller parameters have been 
calculated for sea level and 30,000 ft. altitude, from Appendix I and Ref. 3 
respectively. The results are plotted against equivalent air speed in Figs. 2, 3 
and 4. The very small variation of E, and E> will be observed. 
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FULL ALLOWANCE FOR 
COMPRESSIBILITY 
NO COMPRESSIBILITY EFFECT 
ON ORAG 
——-——- NO COMPRESSIBILITY EFFECT 
ON ORAG OR ON 
PROPELLER EFFICIENCY 


' 30,000 FT. 


250 300 350 400 
Vi M.P.H. 


Fig. 5. 
Example aircraft. Correction to speed at full throttle for temperature variation. 


Finally, the corrections for weight and temperature variations from standard at 
full throttle opening have been calculated from equations (2) and (3). Values of 8 
for use in these equations have been estimated on reasonable assumptions concerning 
the variation of engine power with r.p.m.; the average values are 0.47 at sea level 
and 0.63 at 30,000 ft. The corrections are plotted against equivalent air speed in 
Figs. 5 and 6. 
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FULL ALLOWANCE FOR 
COMPRESSIBILITY 


—-—-— NO COMPRESSIBILITY EFFECT 
ON ORAG. 
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ON ORAG OR ON 
PROPELLER EFFICIENCY 
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Fig. 6. 
Example aircraft—correction to speed at full throttle for weight variation. 


In these figures, the curves which would be obtained if the effects of com- 
pressibility on the drag and propeller efficiency were neglected are also shown. As 
will be seen, the effect of compressibility is to reduce the corrections considerably. 


Figures 7 and 8 show the errors in the corrections produced by errors of +0.1 
in estimation of the various engine and propeller parameters. An error of this 
magnitude in Ep is of no importance, the resulting error in the correction being less 
than 2.5 per cent. of the correction at all conditions. At high speeds, the effect of 
error in II, on the temperature correction is by far the greatest (it has of course no 
effect on the weight correction). The effects of errors in the other parameters vary 
considerably with speed, and become quite large at low speed; fortunately, in this 
region II, and Ey are usually small, and large errors in their estimation are 
unlikely to arise. 
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RECIPROCATING-ENGINE PERFORMANCE REDUCTION 


APPENDIX I 


Formulae for Ql, and Ul, for the conventional four-stroke 
mechanically supercharged engine 


It is known that the performance at constant r.p.m. of the conventional four- 
stroke engine with mechanically-driven supercharger can be represented by 


the relation: 
(P+F) cc (rp, — p/c)/(T,+2). 


Here p, and 7; are the total pressure and temperature in the engine intake, r is the 
supercharger pressure ratio and c the cylinder compression ratio, and z is a 
parameter, dependent on the r.p.m., representing the effective charge temperature 
rise from intake to cylinder, due to compression in the supercharger and heat transfer. 
For p,; and 7, may be substituted 


pi = p(1+4nyM’) with sufficient accuracy, 
T,=t(1+4(y—1)M?) exactly, 
where », represents the intake efficiency and may be taken to be constant. 
Thus 


p[r(1+4nyM?)—1/c] 
t{1+4(y—1) 


(P+F) 


Differentiating at constant p, and remembering that at constant V, and p, M is 
independent of 


Or _ _1+3G-))M (a) 
P+F ot r—1/[c(1+4nyM*)] ar 


It is generally accepted) that the variation of r with temperature at constant engine 
r.p.m. and throttle opening can be represented by 


where C is a constant lying between 0.001 and 0.003. Thus 
roc { 


1+ { 
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Substituting in equation (a) and neglecting 1/c compared with r, 


i, 


where 7; is, as noted already, a function of t and M. This equation applies at 
constant throttle opening: at constant boost pressure, the first term in the square 
brackets disappears. 


To evaluate II, the procedure is as follows. 
At constant pressure, 


M or nyM? (y- 


Therefore 


Hy ~(1+F/P){ )+ 


Ti) l +4nyM? 


at constant throttle opening. At constant boost the formula reduces to 


Il, = —(1+F/P)(y— 1) M?t/(T,+2). 
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AXIALLY-SYMMETRICAL SUPERSONIC FLOW 
NEAR THE CENTRE OF AN EXPANSION 
by 


N. H. JOHANNESEN 
(Mechanical Engineering Laboratory, The Technical University of Denmark) 


and 


R. E. MEYER 
(Department of Mathematics, The University of Manchester) 


SUMMARY 

When a uniform, two-dimensional supersonic flow expands suddenly round a 
corner in a wall it forms a pattern known as a Prandtl-Meyer expansion or centred 
simple wave. If the flow is two-dimensional but not initially uniform, or if it is 
axially-symmetrical, the expansion is still centred, but is not a simple wave. ‘An 
approximate solution is given in this paper for the isentropic, irrotational, steady 
two-dimensional or axially-symmetrical flow of a perfect gas in the neighbourhood 
of the centre of such an expansion. The solution is designed to replace the con- 
ventional method of characteristics in such a region. 


The main application is to a jet issuing from a nozzle that discharges into a 
container with a pressure lower than that in the nozzle; in particular, a formula is 
derived for the initial curvature, at the lip of the nozzle, of the boundary of the jet. 
The solution also applies to the flow near an edge in a boundary wall, and a 
formula is derived for the velocity gradient on the wall immediately downstream of 
the edge. 


1. INTRODUCTION 

The present investigation was stimulated by the research on steam- and air-jet 
ejectors carried out by one of the authors, who was able to show that several 
experimental results concerning ejectors may be explained theoretically if certain 
aerodynamic features of the actuating jet can be calculated”). The experimental 
evidence indicates a correlation between the behaviour of an ejector, with given shape 
of the ducts, and the position of, and the flow conditions at, the point P where the 
boundary of the actuating jet (or rather the mixing zone that develops along it) 
meets the wall of the diffuser (Fig. 1). 


A similar physical problem occurs with a thrust augmenter for aircraft 
propulsion if one of the two mixing gas streams is supersonic. 


Paper received March 1949. 
IThe Aeronautical Quarterly, Vol. Il, August 1950} 
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Fig. 1. 
(1) Actuating nozzle, (2) Diffuser, (3) Jet boundary, 0 Lip of nozzle. 


In order to compute the boundary of the actuating jet by the method of 
characteristics it was found desirable to obtain an approximate solution for the 
field of flow near the lip of an axially-symmetrical, supersonic nozzle. Such a 
solution was therefore sought (Section 2) in the form of an expansion in powers of 
the ratio, R, of the distance from the lip to the distance of the lip from the axis 
(Fig. 2). 


The field of flow in the neighbourhood of the lip of a nozzle may be 
divided into three parts, namely the centred expansion and the regions upstream 
and downstream of it. The first term in a series expansion with respect to R is 
identical with the solution of the corresponding two-dimensional flow problem; to 
this approximation, the flow is uniform in the regions upstream and downstream 
of the expansion, and of the Prandtl-Meyer type in the region of the expansion”. 
The equations for the second term in the expansion are integrated for axially 
symmetrical flow in Sections 3 and 4. From the resulting formule the velocity 
components can be calculated (with an error of the order of R*) and the initial 
curvatures of the jet boundary and of the Mach lines in the expansion region can 
be determined (Section 6). 


The solution applies also to the flow in the neighbourhood of an edge in the 
surface of a body of revolution, or in the wall of a duct of circular cross-section. 
In such cases the velocity gradient on the wall just downstream of the edge may 
be determined (Section 6). 


The corresponding formule for two-dimensional flow are given in Section 8 
(equations (33) to (35)). 


In a jet issuing from an axially-symmetrical nozzle a limit line, and hence a 
shockwave, may be encountered that starts on the downstream border of the region 
of centred expansion. Knowledge of the initial curvature of the jet boundary 
permits an estimate to be obtained for the position at which the shockwave starts 
(Section 7 and Appendix). 
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The calculations are carried through for the case of a perfect gas without 
viscosity or heat conductivity, and the flow is assumed to be isentropic, irrotational, 
steady, and shock-free. The results, therefore, apply only to the case where an 
expansion does in fact take place at the lip of the nozzle or the edge on the surface. 
By introducing the shockwave equations instead of the solution for a centred 
expansion, the method could be adapted to an approximate treatment of the flow 
in the neighbourhood of a lip, or edge, with an attached shockwave. 


Notation 
distance from the axis of symmetry 


ratio of the distance from the centre of the expansion 
to the distance of the centre from the axis 


polar co-ordinates in the meridian plane with origin at 
the centre of the expansion (Fig. 2) 


velocity components in the directions of R and ¢ 
increasing, respectively (Fig. 2) 


velocity magnitude 


angle between the downstream direction on the axis and 
the local stream direction 


local speed of sound 
a, critical speed of sound 
# local Mach angle (sin »=a/q) 


Uy, Vo,@, first terms in the expansions of u, v, a, respectively, in 
powers of R (values at the centre of the expansion) 


u,,v,,@, coefficients of R in the expansions of u, v, a, respectively, 
in powers of R 


q..w velocity magnitude at the centre of the expansion in 
regions of type A and B, respectively 


M = Mach number at centre of expansion (M=q,/a,) 
. . . integrals occurring in equation (14) 
P(u) Busemann’s “ Pressure number ” 
hs focusing parameter of Ref. 4 
A,,C,,C,, K,. 8,, 8, constants of integration 
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2 (a,*/w) cos B, 

rectangular co-ordinate 

adiabatic exponent 

angle of deflection of the streamline at the centre of the 
expansion 

angle at the centre containing the region of centred 
expansion 

A 


limiting value as the centre of the expansion is 
approached along the boundary streamline 


limiting value as the centre is approached along the 
boundary streamline from upstream 


limiting value as the centre is approached along the 
boundary streamline from downstream 


limiting value as the centre is approached along the 
leading Mach line of the region of centred expansion 


limiting value as the centre is approached along the 
final Mach line of the region of centred expansion 


2. EQUATIONS OF MOTION 


Let the radius of the nozzle mouth be taken as the unit of length. We 
introduce a system of polar co-ordinates R, ¢ in a meridian plane, with origin at 
the lip and with measured in the counter-clockwise sense from the normal to the 
axis through the lip (Fig. 2). All the dependent variables are assumed expanded in 
powers of R, with coefficients depending on ¢. Thus, if u and v denote the velocity 
components in the directions of R and ¢ increasing, respectively, 


u(R, ¢)=u, ()+Ru, (9)+. . . 
v(R, ¢)=2, Rv, 


The absence of vorticity is expressed by the equation 


S 
8 
A 
A 
Suffixes : 
| 
| | 
| 
_ 12 (Rv) _ 
Roo R OR =O (2) 
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Fig. 2. 
The equations of continuity and of motion may be combined to give the equation* 


Ris Roo sin ¢—ucos 9)=0. 


R T=. 
(3) 


Bernoulli’s equation is 


9 


where a is the local speed of sound, a, the critical speed of sound, and 


If the series (1) be substituted into (2) and coefficients of the various powers of 
R are equated the following are obtained: 


v,=4u,’, 


etc., where dashes denote differentiation with respect to ¢. Similarly, from (4), 
(7) 


(8) 
. It can 1 be obtained, for example, from the equations (7*), (8) and (10) for motion in general 
orthogonal co-ordinates in a meridian plane in Ref. 3. 
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and then from (3), when use is made of equations (5) to (8), 
(u, + (v,? a,”7)=0, e ‘ (9) 


(v,? + (Uo t+ -D Ug, +(y + 1) + + (Uy? — )u, 


—a,? (v, —u, cos ¢)=0. ; (10) 


Equation (9) shows that the system has two solutions. One, solution A, 
approaches a simple wave with centre at the origin (“ Prandtl-Meyer Expansion,” 
Ref. 2), as R —> 0; the other, solution B, approaches a field of uniform flow. 


3. THE CENTRED EXPANSION (SOLUTION 4A) 
If 


v,7—-a,?=0 . 
(5) and (7) integrate to 


uy (a,/A) sin [A (9+ 8,)], 
(?)=a, cos [A B,)], 


Qo” + V%, =a2(1 [A ,)1) 
(Ref. 2), where 8, is a constant. By (11), equation (10) reduces to 


2u,u,’+ 


and after substituting from (12) this integrates to 


u, 5; (cos x 
{ 1, cos 8, + 1, (n) sin By Al; (n)cos +A, (n) sin +C, } , 


where »=A (¢+8,), C, is an arbitrary constant, and 


1, ()= cos (y/A) (sin y)-* (cos dy, 


I, (n)= sin (y/A) (sin (cos 1/02") dy, 
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1; (n)= sin (y/A) (sin y)~4 (cos dy, 


I, (n)= cos (y/A) (sin (cos i(y-1) dy, 


(cf. Section 6, Table I). The introduction of the constant C, allows any convenient 
value to be chosen for the lower limit, ,, in the integrals. v, is found from (6) 
and (13). 


4. SOLUTION B 


If 
Uy + Vo = 0, 


it follows from (5) that 
u,=w cos 
v,=—wsiny, 
w? =u," + v,? =constant, 
By (15), (6), and (17), equation (10) reduces to 
4 (v,? + + (w? —2a,”)u, —a,? (v, sin @—u, cos . 
and upon substituting (16) we obtain the equation 
(sin? — a,?/w?) u,” —2 cos u,’ +2 (1 —2a,?/w?) u, +S=0, 
where S=2(a,?/w) cos 8,. The general solution of (19) is 
U, (9)=A, cos (24+ 8.)+K,, 
provided the constants A,, K,, 8, satisfy the relation 
A, cos 8, + K, (1 —2a,?/w?)+4S=0, . 
as can be verified by substituting into (19). By (6), 


Vv, (¢)= -A, sin + 
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5. BOUNDARY CONDITIONS 
It will be convenient for the application of these results in the following sections 
to consider as separate branches the upstream and downstream parts of the stream. 
line through the origin. On either of them 9 and @ (the angle that the stream 
direction makes with the axis) may be regarded as functions of R, and near the 
origin, 
@=¢3+(do/dR)R+O (R’), 
6=6, + [(06/0¢) (do/dR)+00/0R] R+O (R*), 
where the derivatives are taken for ¢=¢3, R=0. From Fig. 2, 


and if the streamline belongs to a region where the solution is of the type B, 
00/0¢=0 by (16) and (17); moreover, v (0, ¢3)=v, (@3)=0, so that the curvature of 
the streamline at the origin becomes 

(d6/dR)n-. Uy (On) 
(except for the sign in the case of the upstream branch). 


Similarly, 


q dq/dR=q (éq/00)(do/dR)+q 0q/éR 


on either branch of the streamline; if it belongs to a B region, 0q/do=0, and 


qdq/dR=u, (¢3)+O(R). . . (26) 


It may be noted that v, must be a discontinuous function of # on the border 
between an A region and a B region, by (6), (13) and (18). The Mach line that 
separates the two regions must carry a first order disturbance). However, the total 
derivatives of g and @ along the Mach line must be continuous, and it follows from 
the continuity of either that the limiting value of v, in the B region must exceed 
that in the A region by (q cos »)/(y+1) times the initial curvature of the border 
Mach line, where g and » are taken at the origin in the B region. 


6. INITIAL CURVATURE OF A JET BOUNDARY 


The centre of the expansion is at the lip of the nozzle. Let ¢., 1, ¢r, and 9; 
denote the limiting values of ¢ as R—> 0 on the nozzle wall, on the leading and 
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Fig. 3. 


final Mach lines of the expansion region, and on the jet boundary, respectively 
(Fig. 3). Then 


pit der 
| 
| 
A+ pp =84+ wy, 


where »; and up are the Mach angles corresponding to the ultimate speed, gi, on 
the nozzle wall and to the speed, gp, on the free boundary of the jet, respectively. 
A is the angle at the lip containing the expansion region, 6 the angle of deflection 
of the boundary streamline. 


For 9 <% < 9, the solution is of type B; it is not influenced by the fact that 
the nozzle ends at the origin. Since u,(¢.)=—q1, v) (¢.)=0, we may put w=q,, 
¥=—4,, in this region, whence 


Uy i COS Vo (¢.)= sin (28) 
by (16), (17), and (27). By (20) to (22) and (25) to (27), 
U, (¢,)=2 Cos? py, U, SiN V, (.)+4S 


= dt sin [(dé/dR)r-. - (i cot sin 6, sin? Pu (29) 
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where —(dq/dR)p-,, and —(d@/dR)z-, are respectively the stream direction and 
the velocity gradient on, and the curvature of, the nozzle wall at the lip; these we 
assume to be given. 


For ¢, <¢ < ¢ the solution is of type A; the value of ¢ depends on the 
outside pressure. By (12), (16), and (27), 


(A! tan 
and Bernoulli’s equation may be written 


(gr/ qu)? =(y—14+2 sin? w,)/(y—1+2sin? yy). (31) 


For ¢y < ¢ <¢; the solution is again of type B. Since u, (¢;)=qr, v, (¢;)=0 
we have now 8, =¢;, w=qr, and hence, by (20) to (22) and (25) to (27), 


(d0/dR)z-. —(1/qe) cotur (dq/dR)x-, = — [gr sin 9; sin? vp + (32) 


on the downstream branch of the boundary streamline; 6;=9,;—4* is the stream 
direction just downstream of the centre of the expansion. In the case of a jet, 
dq/dR=0 on its boundary, and (32) gives the initial curvature of the jet boundary. 
If the centred expansion arises from an edge on the wall of an axially-symmetrical 
duct, d6/dR is given on the wall, and from (32) we may calculate the velocity 
gradient on the wall just downstream of the edge. 


For the calculation of the initial slope and curvature of the jet boundary, the 
following steps are necessary. First, u, and u,(¢,) must be calculated from the 
known conditions on the nozzle wall, immediately upstream of the lip (equation (29)). 
Secondly, gp must be found from the known pressure outside the jet and the initial 


Fig. 4. 
(1) Leading Mach line, (2) Final Mach line, (3) Jet boundary. 


* A table of P (u)=1,090° —u—(180°/7) A-! cot—! (A-! tan is given in Ref. 5. 
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TABLE I 


sin I, (n) I, (n) (n) I,(n) 
.73167 47°.026 — .025191 056596 053642 — .023868 
.74348 48.029 —.054656  .115879 .107910 — .050833 
.75530 49.052 —.089006  .178076 162857 081169 
76712 50.096 — 128971 243441 .218536 115202 
.77894 51.163 —.175440  .312256 .274999 153318 
.79076 52.256 —.229498  .384842 332301 195981 
80257 53.377 —.292490 461560 .390495 — .243748 
81439 54.527 —.366095  .542815 449630 — .297298 
82621 55.711 —.452442  .629062 509751 — .357467 
83803 56.933 —.554272  .720812 570885 — .425293 

for sin n.=.71985 (n,=46°.0420) and y=7/5 (for air). 


state of the gas. sup and 4 can then be calculated from (31) and (30), and 4 follows 
from (27).* wu, (?e) is found from (14), where the integrals have to be evaluated 
between the limits »,=cot~'(A~'tanp,) and 7=,+AA. Finally, the initial 
curvature of the jet boundary is given by (32), with (dq/dR)k-,=0. 


The integrals that occur in (14) have been tabulated in Ref. 1 for 
46°.0420 = n <= 56°.9326—these values correspond, for air with y=7/5, to the 
Mach numbers 2.7300 and 3.8928—and the table is reproduced in Table I; it covers 
all cases of expansions with initial and final Mach numbers between those limits. 


7. THE POSITION OF THE POINT AT WHICH THE SHOCKWAVE 
STARTS 


In Ref. 1 the computation of the field of flow has been carried further in a 
particular case by the method of characteristics. The nozzle was taken to generate 
a pure source flow in a cone of semi-vertex angle 4°, with a Mach number 2.7300 
at the lip (so that ¢,= —86° and ¢,=72°.51; Fig. 4). A centred expansion to 
M=3.893 is then. found to lead to a deflection of the boundary stream- 
line through 6=20°.074 (A=26°.68 and ¢,=99°.19; it was also found that 
and No use. was made of 


* In the case of an edge in the wall, 5 is given and yr is found by the help of the P, »—table, for 
P (ut)—P (ur) by (27) and (30)(?). 
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equation (32) in computing the boundary of the jet; this equation predicts an initial 
radius of curvature | dR/dé@|=10.175, a value which is found to agree well with 
the results of the computation of Ref. 1. 


It is of interest to examine how far the results of the present. investigation, 
together with those of earlier publications, permit an estimation of the position at 
which the first shockwave will start in a jet of this type. The considerations will 
be based on the well-known conjecture“) that the shockwave starts either where a 
limit line starts, or at the cusp of a limit line. It is to be expected that the region 
of the centred expansion will not contain a limit line. If the downstream border 
of this region, i.e. the final Mach line, meets a limit line the position of its limit 
point can be determined by a method indicated in Ref. 4, where a length parameter 
hg is introduced that vanishes at the limit point and a differentia] equation is given 
that governs the variation of hs on a Mach line (equation (43) of Ref. 4). To 
integrate this “ focusing ” equation it is necessary to know (i) the initial value of hy 
on the downstream side of the final Mach line, at the centre of the expansion, and 
(ii) the velocity distribution on this Mach line, on which the coefficients of the 
equation depend. The former is given exactly by the present theory, but the velocity 
distribution remains unknown except for the values of the velocity components and 
their derivatives at the lip. In the absence of further information estimates can only 
be substituted for the coefficients in the focusing equation of Ref. 4, based on their 
known values at the lip. 


In this way, two estimates are derived in the Appendix for the position of the 
point where the shockwave starts. For the special case of Ref. 1, equation (Al) 
predicts a limit line starting where r=1.35 on the final Mach line (r denotes the 
distance from the axis). The more refined estimate (A4) gives r=1.315. In Ref. 1 
the first intersection of adjacent Mach lines of the characteristic network is found 
to occur at Q (Fig. 4), where r=1.426 on the final Mach line of the expansion 
region. This result, however, depends on the step size employed for the numerical 
process, and from it it may be concluded only that the limit line starts on this Mach 
line somewhere between the points where r=1* and r=1.426. 


8. COMPARISON WITH THE TWO-DIMENSIONAL CASE 


For the two-dimensional flow near the centre of an expansion the equations (1) 
to (13), and (15) to (28), remain valid if the last term on the left hand side is deleted 
in equations (3), (10), (13), (18), (19), and (21). Hence, the solution B in the two- 
dimensional case is obtained from that in Section 4 by putting S equal to zero, and 


* The limit line cannot start at the lip itself. For if the lip were a point of an envelope of the 
Mach lines of the region downstream of the centred expansion, the initial curvature of the 
jet boundary would be infinite, and this is not the case. 
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instead of equations (29) and (32) we have 
and (d6/dR)n=o — (1/ qr) cot (dq/dR)p-o= — u, (¢r)/ (Gr Sin 
respectively. For solution A we now find 
u, ()=C, [sin (cos YP, 


instead of (14), where C, is a constant of integration. In two-dimensional flow the 
expression 


p=0- | (1/q)cotu dq 


is a characteristic parameter'*), and equations (33) to (35) show that its derivative 
in the stream direction increases from the upstream side to the downstream side of 
the centre of an expansion in the same ratio as the expression 


[sin » (cos }/(q sin 2x), 


where »=cot~' tan 


The equations (32) and (34) differ by a term which is proportional to the radial 
velocity component just downstream of the lip; if this component of the velocity 
is positive the additional term reduces the initial curvature of the jet boundary, and 
hence encourages the formation of a shockwave in the case of axially-symmetrical 
flow (cf. Appendix). On the other hand, a low value of u,(¢r) is seen to discourage 
the breakdown of potential flow. Equation (13) shows that an increase (decrease) 
in the magnitude of u, in the transition through the centred expansion is common to 
two-dimensional and axially-symmetrical flow if M<(>)¥V { 4y/(y-—1)}; if the 
radial velocity is positive, the last term in (13) favours an increase in u,, and thus 
hastens the formation of a shock. A positive radial velocity at the lip just upstream 
of the expansion, however, is seen to reduce the initial value of u, in the expansion 
region (equation (29)). 
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*2cos times this derivative of 8 is increased by the reciprocal of the growth factor which 
governs the steepening of wave fronts in two-dimensional flow (°). 
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APPENDIX 
(The notation used is that of Ref. 4) 


Estimates for the position of the point where the shockwave starts 


Additional symbols not included in the main notation given on pages 129 and 130 
e 6- 
Rs, codins of curvature of the free boundary at the lip 
r logr 
y2 g 
n? g-e 
ks curvature of Mach line 
If the curvature of the final Mach line of the expansion region and the variation 
of the velocity components along it are neglected, i.e. if equations (44) of Ref. 4 are 
assumed to hold, equation (51) of Ref. 4 remains valid; the streamline curvature at 
the lip, as defined in Ref. 4, is —(d0/dR)g-,. If, in addition, the radial velocity 
component on the final Mach line is neglected, so that (44a) holds, but it is not 
assumed that this Mach line is inclined to the axis at the Mach angle, then 
equation (52) of Ref. 4 remains valid, provided N is replaced by 


N’= —(y+1)/ { cos?u sin » sin (@—»)}, 


where @— =e is the inclination of the final Mach line. Thus hy vanishes, and the 
shockwave may be assumed to start, where 


on the final Mach line, with 
Rs, (dR/d®)p-,. 

If it is assumed that equations (44) but not (44a) of Ref. 4 hold on the final 
Mach line, i.e, that the radial velocity component is constant on it and equal to its 
value at the lip*, equation (83) of Ref. 4¢ reduces tot 

d’h 4 2e+1 dhs 
dr* 


& 


* It should be noted that neither approximation is consistent in the sense that a characteristic 
with the properties assumed could exist. For it is shown (7) that if a Mach line in isentropic, 
irrotational, axially-symmetrical, steady, supersonic flow is straight and if the velocity 
components are constant along it, then the radial velocity component must vanish on it, and 
its inclination to the axis must be equal to the Mach angle. 


t In that equation, the first term in the square bracket should read 27,R/v,"; and the left hand 
side of equation (84) should read 9R/hsd8. The definition of A given in Ref. 4 differs from 
the one adopted here. 


t Since a single Mach line is being considered, it can be assumed without loss of generality that 
0h, /da=0. 
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sin 6 [a5 2 


with e= — 


sin? 6 


sin 
8” 4 cos? sin? 


COS 
(A2) can be written 

d*y/dr* +(g—e*) y=0, 
with y=rhs and r=logr. If v’=e?—g>0O*, hy vanishes for 


dy/dr—vy 


If n?=g—e?>O0f, hs vanishes for 


nlogr= —tan 


by (3) of Ref. 3, and it may be assumed without loss of generality that hs=1 (and 
hence y=1) for r=1 on the Mach line. Since it is a “plus Mach line,” the upper 
sign must be taken in equation (51) of Ref. 4 and «=@—j, whence 


(dy/dr),-, =e-—(y+ 1)/ { 2 cos? psin sin Rs, } 
where 1/Rs,=—(d0/dR)z-, 


in the notation employed in Sections 5-7 above. 


* In the case treated in Ref. 4, g=0 and e= —1/4. 
+ This sign is found to hold in the case discussed in Section 7 above. 
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NOTE ON THE DEPENDENCE OF FLAP HINGE 
MOMENT DERIVATIVES ON HINGE POSITION 


by 


W. J. DUNCAN, D.Sc., F.R.S., F.R.Ae.S. 
(Professor of Aerodynamics, College of Aeronautics, Cranfield) 


SUMMARY 

It is shown to be a consequence of kinematics and statics that the aerodynamic 
derivative expressing the rate of change of hinge moment of a plain flap with the 
angular setting of the flap varies parabolically with the chordwise position of the 
hinge. The derivative expressing the rate of change of hinge moment with incidence 
of the main surface varies linearly with the hinge position. Similar conclusions 
hold for the derivatives of the hinge moment with respect to the corresponding 
angular velocities. 


The object of this note is to obtain all the information on hinge moment 
derivatives which can be obtained from kinematics and statics alone. 


Take the main surface and flap in a certain datum configuration, placed in a 
uniform airstream of fixed characteristics and let the hinge be at a datum position 
in the chord line. For a given incidence of the main surface the aerodynamic forces 
on the flap will be statically equivalent to a force N, normal to the chord and 
through the hinge axis, a tangential or chordwise force T and a hinge moment H. 
The force T is small and will at present be neglected but its influence on the hinge 
‘moment would not be negligible if the hinge axis had a large vertical offset from 
the chord line. WN is taken positive downwards and H positive in the trailing edge 
down sense. In general, N and H will depend on the angular deflection € and on z, 
the normal downward displacement of the flap*, and we shall express them by 
means of derivatives as follows 


(1) 


(2) 


Now, with the flap in its datum configuration, let us shift the hinge axis rearward 
along the flap chord through the distance /, and then let the flap be deflected through 


*The reason for introducing this will become apparent later. 
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the angle ¢. When ¢ is small this raises the original hinge by the amount [¢, 
ie. z= —I€. Hence N and H, still referred to the original hinge axis, are given by 


aN 
N=N,+é 


By statics the normal force N’ and hinge moment H’ referred to the new hinge axis 
are given by 


(3) 


(5) 


Hence by (3) and (4) 


(8) 


_ aH ON) , 


Thus the hinge moment derivative varies parabolically with hinge position. This 
simple fact will be obscured when the conventional non-dimensional coefficient b, 
(the rate of change of hinge moment coefficient with £) is calculated, on account 
of the division by the area of the flap aft of the hinge, since this varies as the 
hinge is shifted. 

Suppose next that A is some parameter affecting N and H but not involving 


change of position of flap relative to main surface. Then we derive immediately 
from (5) and (6) that 


(9) 


and 


| 
| 
ON’ ON F 
aN’ aN 
0H’ 0H 
In particular, if we identify 4 with the angle of incidence 2 we derive ; 
OH’ AN 
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Hence this derivative varies linearly with hinge position. Again this simple result 
is masked in the conventional non-dimensional coefficient b, (the rate of change of 
hinge moment coefficient with 2). Similar arguments apply to the damping 
derivatives. Thus, corresponding to (8) we have 


(12) 


0H’ 0H (2H _ aN 
dg =a az ) 


while (11) corresponds to 


All the results given apply to flaps in two and three dimensions. 


The author is not aware that the foregoing simple arguments have been 
published previously in an explicit form. However, the nature of the dependence 
of the derivatives on the chordwise position of the hinge is implicit in the work 
of Kiissner and Schwarz") and of W. P. Jones‘) on flutter derivatives. 


It is suggested that it may be a convenient technique in wind tunnel work on 
hinge moments to measure H and N and their dependence on € and z for a single 
chordwise hinge position rather than to use a succession of chordwise hinge positions. 
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ON THE CHORDWISE LIFT.DISTRIBUTION AT THE 
CENTRE OF SWEPT WINGS 
by 


D. KUCHEMANN 
and 
J. WEBER 
(Aerodynamics Department, Royal Aircraft Establishment) 


SUMMARY 

The present note deals with the peculiarities of the lift distribution on swept 
wings. The influence of the angle of sweep is expected to be greatest in its effect 
on the chordwise lift distribution at the centre section of the wing. It is assumed 
that the essential characteristics of this chordwise lift distribution are predominantly 
determined by the special configuration of the bound vortices and not so much 
by the trailing vortices. In calculating the downwash, induced by a continuous 
chordwise distribution of swept vortex lines of infinite span, it is found that the 
downwash at the centre line is to be expected to differ from its two-dimensional 
form mainly by a term proportional to the local vortex strength. Solutions of this 
simplified downwash equation have been obtained. They differ from the ordinary 
flat-plate distribution only in the exponent which instead of 4 becomes 4-9/7. It is 
suggested that this type of function is more appropriate for the centre-section of 
swept wings and that it may replace the ordinary flat-plate distribution in all those 
cases where the lifting surface is swept. 


1. INTRODUCTION 

It is known that both the spanwise and the chordwise lift distribution on 
swept wings show appreciable differences from those on straight wings. The 
chordwise lift distribution in the centre region in particular differs essentially from 
the distribution on two-dimensional aerofoils. Generally, the chordwise lift farther 
out on the wing resembles that on unswept two-dimensional wings and, for 
uncambered aerofoil sections, the “flat-plate” distribution is known to give a 
good first approximation. At the centre-section, of swept-back wings, however, 
the local lift in the front part of the section, as shown by experiments, is reduced 
compared with that of the flat plate; whereas the rear part shows a somewhat 
higher local lift. This generally results in a reduction of the local lift coefficient 
and also in a rearward displacement of the aerodynamic centre. These changes 
are greatest in the plane of symmetry. The flat-plate distribution cannot give a 
sufficiently accurate first approximation in this case and more terms of the Birnbaum- 
series for the chordwise lift distribution are needed if the ordinary method is to 
be applied. Attempts should be made therefore to find a function whose first 
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term would give a satisfactory first approximation. The present note suggests a 
type of distribution which might be appropriate and which may thus replace the 
flat-plate distribution in the plane of symmetry in all those cases where the lifting 
surface is swept. 

In order to simplify the mathematical approach it is assumed that the essential 
characteristics of the chordwise lift distribution are predominantly determined by 
the special configuration of the bound vortices and not so much by the trailing 
vortices. The trailing vortices are neglected therefore; and a swept wing of infinite 
span with a constant spanwise lift distribution is taken as a basis. 

To calculate the lift distribution at the centre of the wing the downwash field 
of a given vortex distribution is first needed. In determining the downwash, 
induced by a continuous distribution of swept vortex lines along the chord, it is 
found necessary to avoid the usual simplification of calculating the value on the 
chordline; the proper limit as z (the vertical co-ordinate) tends to zero has to be 
taken. Since one of the downwash terms tends to infinity at the chordline, a further 
simplification is suggested tentatively, allowing either for the finite thickness of the 
wing profile or for a suitable rounding off of the vortex lines at the centre. The 
downwash is assumed to differ from its two-dimensional form only by a term 
proportional to the local vortex strength. Solutions to the simplified downwash 
equation are then obtained (equation (13)). Values for one of the parameters in 
this solution are determined by using some experimental data. Finally, an expression 
for the chordwise lift distribution at the centre is obtained, where only a general 
factor, i.e. the total lift coefficient, remains to be determined. This must be left to 
a more complete investigation where the effect of the trailing vortices is included. 
But the type of distribution found agrees fairly well with experimental data and 
appears to cater for the effect of the angle of sweep. 


Notation 
rectangular co-ordinate system; x in direction of flight; y to starboard; 
z downwards 
wing chord 
distance of aerodynamic centre from leading edge 
angle of incidence 
angle of sweep 
velocity of main flow 
velocity increment in x-direction 
induced velocity in z-direction 
strength of isolated vortex line 
vorticity of continuous vortex distribution 


difference of the pressure coefficient on upper and lower surface of 
the wing 


local lift coefficient 
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parameters, independent of x 

a number between 0 and 1 

numbers, dependent on n (equation (9)) 
variable of integration 

a finite positive number, small compared with 1. 


2. CALCULATION OF THE DOWNWASH 
A standard co-ordinate system has been used: x in the direction of flight, 
y to starboard, and z vertically downwards (Fig. 1). An isolated infinite vortex 


line in the plane z=0 with a kink at x= y=0, swept back at an angle 4, induces a 
downwash velocity v, in the plane y=0, given by 


x [1 ~sin | (1) 
V~, cos x? +(z/cos (x? + 2?) 


I is the strength of the vortex line; V, the velocity of the main flow. The down- 
wash of a continuous distribution of such vortex lines of vorticity y (x) between 
x=-—1 and x= +1 is obtained conveniently by integration in three parts as follows: 


+1 


dx’ 


where I, (x, 0, 


+1 


I, 0, z)= —sin of 


ax, 
[ (x— +27} 


y (x) (z/cos 6)? , 


The evaluation of these integrals requires some special care since the denominators 
of the integrands vanish at the point x’=x, z=0. 


Putting z=0 before integrating is permissible for the first integral 


+1 


dx 


and the principal value has to be taken. This part corresponds to the well known 
expression for the downwash of a two-dimensional vortex distribution. 


In treating the third integral, the range of integration is split into three parts: 
from —1 to x—a; from x-—a@ to x+a; and from x+a to +1; a being a finite 
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Fig. 1. 


positive number, small compared with |. In the first and last of these intervals 
the integrand has no singularity; therefore the limit as z tends to zero can be taken 
before integrating. Consequently these parts do not contribute to the downwash. 
Within the remaining interval -a<x—x’<+a it is assumed that y(x’) is a 
continuous function which can be replaced by its value at one point within the 
interval, at x’=.x say. I, can then be written as 
+a 
1, (x, 0, 0)=sin y (x) im| 


-a 


{ (x—x’P +(z/cos } x’ 


, 


(x — x’), 


which gives 
(1 + (z/a/*}*—sin » 
[1 


I, (x, 0, 0)=y (x) lim log, 


so that finally 
I, (x, 0, O)== f(@).y@) 


where 


This result does not depend on the actual value of a. The double limit z—- 0, 
a—> 0 gives the same result. Obviously z must tend to zero more rapidly than a, 
so that z/a—>0. Otherwise the contributions of the first and third interval could 
not be ignored. 

This additional downwash term depends on the local value of y (x) only, and 
corresponds exactly to the velocity increment v, at the centre line of a thick swept 
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wing at zero incidence (see e.g. F. Ursell)). The additional downwash will spread 
over the centre region of the wing, and will fade out sideways. Its local value will 
depend on the distance from the centre line in terms of the wing chord. This 
additional downwash disappears when ¢=0. 

It is the second term I, in equation (3) which requires special treatment, since 
I, shows a logarithmic infinity for z—>0. Such a singularity cannot occur in 
practice. It is caused by the assumed discontinuity of the vorticity vector at the 
centre. 


3. SIMPLIFICATION OF THE DOWNWASH EQUATION 


One way of getting a reasonable approximation for the second downwash term 
is to calculate the downwash on the surface of the aerofoil of finite thickness. 
This leaves the singularity of the vorticity vector inside the profile section. 
Assuming a special distribution of bound vortices, the induced downwash has been 
calculated along a line above and parallel to the centre line. The result in Fig. 2 
shows that, in this case, both the second and the third term are roughly proportional 
to the vortex distribution chosen. This result appears reasonable since the infinity 
of I, at z=0 indicates that the vorticity in the immediate neighbourhood of the 
pivot point is of paramount influence. 

A similar result will probably be obtained by rounding off the vortices at the 
centre, thus removing the discontinuity of the vorticity vector. Only a closed 
additional vortex of a strength proportional to y(x) of limited extent is needed to 
modify the kinked vortex line into a rounded one. Again, a finite value of the 
downwash will be obtained where one part should be proportional to the local 
vortex strength. An application of the method developed by Ursell’) to this case 
might provide a proper solution to this problem. 

In this note an approximate solution is attempted. Although the above 
estimates are not conclusive, they suggest that a downwash equation of the following 
form may be considered tentatively. 


+1 
2nV 


+ay(x) 7 at y=0, (6) 
1 


where o is a factor, dependent not on x but on the angle of sweep ¢, and probably 
on the thickness ratio and shape of the aerofoil section. It is possible that, for 
very thin sections and kinked vortex lines, the value of o tends to infinity. This 
would link equation (6) to the general equation (2). o=0 represents the case of 
the (two-dimensional) infinite sheared wing. 

It is shown in Section 5 that the solutions of equation (6) agree with 
experiments, thus corroborating the suggested type of downwash equation. 
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= 45° 
y=0 


4. SOLUTION OF THE DOWNWASH EQUATION 


The integral equation (6) has solutions of the following type, 


with y—> 00 for x= +1 (leading edge), and y=0 {or x= — 1 (trailing edge). = is a 
factor, independent of x, and 0<n<1. That y(x) is a solution of equation (6) 
can be seen from the fact that 
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0-4 
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Fig. 3. 


where K, and K, are certain constant numbers dependent only on the 
value of n, and 


The transformations 
l-x 
have been used. Numerical values of K, and K, are plotted in Fig. 3. In particular 
K,=-"= ; K,=0 for n=4 (two-dimensional flat plate), 


K,=-t/2; for n=}, 
; for n=0. 


5. APPROXIMATE VALUES FOR THE EXPONENT n 
It is doubtful whether reasonable values for the parameter + can be 
obtained from equations (6) and (7) since the effect of the trailing vortices has been 
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Fig. 4. 


ignored, so that the total lift coefficient cannot be determined this way. But the 
shape of the chordwise lift distribution, that is the value of n, can be estimated 
more easily by using some experimental data. This estimate makes use of the 
simple fact that the aerodynamic centre of the vortex distribution (7) is given by 


 n=1-2h/c, ‘ . (10) 


where fA is the distance of the aerodynamic centre from the leading edge. 
Experimental results in Fig. 4 suggest that 


is a sufficiently accurate approximation for the relation between the aerodynamic 
centre position at the centre-section and the angle of sweep. 


Combining equations (10) and (11) 
n=4 for the unswept (or infinite sheared) wing (flat plate), 
0<n<4_ for the centre of swept-back wings, 
4<n<1_ for the centre of swept-forward wings. 
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ERRATA 


The function 


y (x)=constant x ( its 


(13) 


represents the type of the chordwise lift distribution at the centre of swept wings 
fairly well, considering that the approximation is of the same order as that of the 
flat-plate distribution in the two-dimensional case. A comparison with experimental 
results in Fig. 5, where a constant factor has been applied to give the right value 
of C,, shows reasonable agreement. 


It is to be expected that, in moving outwards from the centre, the chordwise 
lift distribution will change gradually from that given by equation (13) at the centre 
to a flat-plate distribution, at about one chord from the centre, i.e. from n=4—9/7 
to n=4. The wing tips will show a similar effect in the opposite direction, i.e. 
n> 4. 
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Amendments to “A Note on the Bending Moment Induced in the 
Booms of a Spar at the Point of Apphcation of a Concentrated 
Load.” (Part IV, Volume I, February 1950). 


10th May 1950 


I regret that in my note on the above subject, published in The Aeronautical 
Quarterly, February 1950, some errors have come to light. 


I would be most grateful if the corrections could be given publicity, and I am 
accordingly enclosing a statement of the corrections. 


H. F. WINNY, 
Head of Technical Department, 
The Fairey Aviation Co. Ltd 


Page 285. Line 5. Replace 32?—? by 


and 2? — 3? by 
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Line 18. Replace ( 


a?+ (2? + B?)? 
aB ) by ap 
Line 19. Replace 6?) by (2*+ 8?) 

a? + B? a? + B? 


and K in numerator by K°. 


Line 22. Replace 


Replace existing Fig. 2 by new Fig. 2. 
Line 1. Replace 2/8 by 8/z. 


K+1 K+1 
Line 3. Replace by 4 


Line 3. Replace K in numerator by K*’*. 
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Fig. 2. 
Maximum bending moment carried by boom. 


Page 287. 
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